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Abstract 

CP violation from physics beyond the Standard Model (SM) is investigated 
in the Cabbibo-suppressed decay r — > Kixv T with polarized r leptons, to which 
both the J p = 1" resonance K* and the J p = 0+ resonance Kq contribute. In 
addition to the CP-odd rate asymmetry, r polarization enables us to construct 
three additional CP-odd polarization asymmetries that can be enhanced due to the 
interference between the K* and Kq, and whose magnitudes depend crucially on 
the Kq decay constant, fx*- Taking a QCD sum rule estimate of fx* = 45 MeV 
and the present experimental constraints on the CP-odd parameters into account, 
we estimate quantitatively the maximally-allowed values for the CP-odd rate and 
polarization asymmetries in the multi-Higgs-doublet (MHD) model and the scalar- 
leptoquark (SLQ) models consistent with the SM gauge symmetry where neutrinos 
are massless and left-handed as in the SM. We find that the CP-odd rate and 
polarization asymmetries are of a similar size for highly-polarized r leptons and, 
for their maximally-allowed values, CP violation in the MHD model and two SLQ 
models may be detected with about 10 6 and 10 7 r's at the la level. 
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The decay of the r, the most massive of the known leptons, can serve not only as a 
useful tool in the investigation of some aspects of the SM but also as a powerful exper- 
imental probe of new physics phenomena |T[. One phenomenon where new physics can 
play a crucial role is CP violation. In light of this aspect the r decays into hadrons have 
been recently studied as probes of CP violation in the scalar sector of physics beyond the 
SM@|1|. 

In the present paper we extend the previous work |5j by Choi, Hagiwara and Tanabashi 
to probe CP violation in the Cabbibo-suppressed r decay r — ► Kitv T with polarized r's. 
The decay mode is dominated by the contributions of the two lowest vector and scalar 
resonances, K* and Kq, with different spins and relatively large width-to-mass ratios |J: 

K* : J p = l", m K *= 892 MeV, T K * = 50 MeV , ( ) 

K*: J p = + , m*. = 1430 MeV, r K * = 287 MeV, W 

and the mode is expected to have larger scalar contributions than the 2ir or 3ir modes 
due to the s quark mass much larger than the d quark mass. In the light of these aspects, 
the Cabbibo-suppressed mode is worthwhile to be investigated in detail. 

Including possible contributions from new physics with massless left-handed neutrinos, 
we can write the matrix element for the decay r~ — > (K-k)~v t in the general form 

M = V2G F [(1 + x)u(K -)-fP-u(p, a) J M + r)u(k, -)P+u(p, cr)J s ] , (2) 

where P± = (l±7s)/2, Gp is the Fermi constant, and p and k are the four momenta of the 
r lepton and the tau neutrino, respectively. The parameters x an d tj, which parametrize 
contributions from physics beyond the SM, are in general complex. In eq. (0), the helicity 
of the t~ is denoted by o (o = ±1) with its spin quantization direction along its neutrino 
momentum direction. The tau neutrino is assumed to be massless and left-handed as in 
the SM so that the helicity value is —1/2 as indicated by the negative sign in its spinor 
u(k, —). The vector and scalar hadronic matrix elements 

J M = sin 9 C ({Kit)~ | s-f^u \ 0), 

J s = sin 6 c ((Kn)-\7su\0), (3) 

with sin^c = 0.23 for the Cabbibo angle 8q are related through the Dirac equations to 
the s and u quarks at the quark level and their explicit form can be parametrized in terms 
of two form factors Fk{q 2 ) and Fs{q 2 ): 



F K (q 2 ) ~ qj jA (gi - ^ + ^C K F s (q 2 )q, 



Jf, = \/2sm6c 

J s = V2sm6 c ( ^ ) C K F s (q 2 ), (4) 
1 m e — m„ 



where q\ and q^ are the four-momenta of tc and K, respectively, m s and m u the s and u 
current quark masses, and q is the four- momentum of the Kti system; q = q 1 + q 2 . 
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The coupling strength Cx denoting the scalar contributions is determined by the 
Kq decay constant fx* and the coupling strength Qk*k-k of the Kq to Kit under the 
assumption of Br(A" * -A Kn) = 100%: 

C K = fK f K « K \ (5) 



v3m 

The value of gx*K-w is 4.9 GeV from the measured Kq — > fT7r decay width IX/Tq — > Kn) ~ 
287 MeV. Even though the Kq decay constant fx* is not experimentally measured, it 
has been estimated by several model-dependent methods. The QCD sum-rule estimate 
in the Kq narrow-width approximation is fx* ~ 31 MeV @ and an effective Lagrangian 
estimate of fx* including the width effects leads to a larger value of about 45 MeV ||, 
which approaches the pole dominance result of 50 MeV |§. In light of these present crude 
estimates of the Kq decay constant, we adopt for the Kq decay constant the effective 
Lagrangian estimate |§ 

f K * = 45 MeV, (6) 
anticipating that the decay constant fx* is measured more precisely in future experiments. 
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Figure 1: Definition of the angular variables in the r and (Kit) rest frames for the Cabbibo- 
suppressed r decay t~ — > (Kn)^u T . The (JOr) momentum direction in the r rest frame is 
denoted by the polar angle G and the azimuthal angle $ with respect to the r momentum 
direction defined as the positive z-axis, and the K momentum direction in the (Kit) rest 
frame denoted by the angles 9 and (f> with respect to the same positive z direction. 



In spite of the several resonance contributions to Fx and Fg, the form factors are 
approximated to be the propagators of the lowest-level resonance states K* and Kq\ 

F K (q 2 ) = B K *(q 2 ), F s (q 2 ) = B K *(q 2 ), (7) 
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respectively, where Bk* and Bx* are parametrized in the Breit-Wigner form with the 
momentum-dependent widths fT^f : 



B 



rn 



x \ 



x 



m x — Q 2 — irrixT 



x \ 



T x (q 2 ) = T x D x (q 2 



(8) 



for X = K* or Kq. We adopt for the momentum-dependent widths the parametrizations 



of the r-decay program package TAUOLA JTO] 

1 3 

! m K * Pk(Q z ) 

D K *(q 2 







for q 2 > {rriK + m^) 2 
for q 2 < (mx 1 "' ,J 



m K » 



D f 







for q 2 > {m K + m w ) 2 
for q 2 < (m K + m v ) 2 



(9) 
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Let us now calculate the helicity amplitudes of the Cabbibo-suppressed r decay. In 
general, the r three-body decay into (K-k)v t is described by five independent kinematic 
variables. Because the hadronic (Kir) system is solely determined by the lepton momen- 
tum transfer, it is convenient to consider two reference frames, namely the r rest frame 
and the (Kir) rest frame as shown in Fig. [l|. We define the momentum direction of the 
virtual K* or Kq in the r rest frame by the polar angle G and the azimuthal angle $ 
with respect to the r momentum direction taken to be along the positive z-axis and the 
K momentum direction in the (Kir) rest frame by the angles 9 and 0. The rotational 
invariance of the total system with respect to the r momentum direction allows us to 
take $ to be zero in the calculation of the helicity amplitudes, while the azimuthal angle 
$ is employed to describe r polarization, especially the transverse polarization of the r. 
Setting $ to be zero and using the 2-component spinor technique , one can obtain the 
helicity amplitudes M„ of the Cabbibo-suppressed decay r~ — > (Kt[)~v t in the reference 
frames defined in Fig. |l] as 



M+ 



2G F sm9 c (l + x) 



m 2 _ g2 



m T cos — 



-2P K \ cos 9 cos — + sin 9 sin — e irp + 
I 2 2 



M_ 



2G F sin# c (l + x) 



mi 



m T sin 



+2Pk I cos 9 sin — — sin 9 cos — e 
I 2 2 




+ C )C K F s (q 2 



\ o ©1 




1 COS 17 COS — 


> F K (q 2 ) 





+ i)C K F s (q 2 



- 1 



9 • 9 
cos v sin — 
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where d is the angle between the K momentum in the (Kir) rest frame and the (Kir) 
momentum in the r-rest frame, i.e. cosi? = sin 9 cos <fi sin O + cos 9 cos 0, and Pk is the 
size of the K momentum in the virtual K* and K£ rest frame; 

P K (q 2 ) = ^X 1/2 (q 2 ,rnl,m 2 K ), (12) 
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with X(x, y, z) = x 2 +y 2 + z 2 — 2xy — 2yz — 2xz. The dimensionless parameter £ determines 
the relative size of the scalar contributions to the vector ones: 



m 



K* 



(m s -m u )m T \1 + x J ' 



(13) 



Clearly, the initial r~ system is not CP self-conjugate and so any genuine CP-odd 
observable can be constructed only by considering both the Cabbibo-suppressed r~ decay 
and its charge-conjugated r + decay, and by identifying the CP relations of their kinematic 
distributions. Before constructing all the possible CP-odd asymmetries explicitly, we 
calculate the helicity amplitudes for the charge-conjugated process r + — > (i^7r) + z/ r . The 
helicity amplitudes M 5 of the r + decay in the same reference frame as in the r~ decay 
are given by 



2G F sin# c (l + x*; 



mi 



6 /mi, 
m.sin-j-^ + rl^F^g 2 ) 



e 



e 



-2P K {-cos 9 sin h sin 9 cos — + 



Q 



1 cos i? sin 
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M_ 



-2G F sin0 c (l + x* 



mi 



m T cos ■ 



6 (m 2 K , 



F K (q 2 



Q 



C C K F s (q 2 



[ - 6 © - / m 

-2P K \ cos 9 cos — + sin 9 sin — e~ %4> + 

I 2 2 \Vr 



— l^j cos $ cos -^-| 



.(14) 



It is easily shown that, if the parameters 7] and \ are real, the helicity amplitudes (|TT|) 
for the r~ decay and (0) for the r + decay satisfy the CP relation: 

M±(8; g 2 ; 0, 0) = TM T (9; g 2 ; 0, -0). (15) 

With the results in eqs. (O) and (TJ) for the r 1 * 1 decay helicity amplitudes, one can 
describe the decay of an arbitrary polarized r lepton by superposing the two helicity 
states. Generally, a pure r polarization state, which is polarized in the direction (9 p ,(f) p ), 
is given in terms of the helicity states by 



10, 







pi rp 



) = cos-jf |+) + sin-jfe l * p |-) 



(16) 



where P and 4> p are the polar and azimuthal angles in the r rest frame, respectively, with 
respect to the r momentum direction in the laboratory frame. In this case the decay 
amplitude to the final state |0, $) of the virtual K*/Kq system is expressed in terms of 
the helicity amplitudes as 

(e,$|0 p ,0 p ) = (e,o|0 p ,0 p 



$) = cos + sin ^V ( ^- $) M_, 



(17) 



where the rotational invariance of the total system with respect to the r momentum 
direction has been used to relate the first and second expressions. It is noteworthy that 
only the difference $ — <p p between the azimuthal angles $ and 4> p appears. 
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Taking into account the polarization degree P for a partially-polarized r beam, one 
can decompose the differential decay rate into four independent terms: 

dv = ~d(r ++ + r__) + ~Pcos0 p d(r ++ -r__) 

+Psinfl p cosOp - $)Re(cffV) - Psin# p sin(0 p - $)Im(dr+_), (18) 
where the helicity-dependent terms are defined as 

dT aa , = - \ r 1 ( 1 - -^r ) M a M*,P K d x ff dcos6 deb dcosQ (19) 
(27r) J 32m T y m l T J v 

For the sake of notational convenience, we introduce d<& 3 for the phase space element 
d^fq 1 d cos 9 deb d cos and denote the four independent terms in eq. (|ll|) as 

dVi rf(r ++ + r__) dr 2 d(r ++ -r__) 



rf$ 3 d§ 3 rf$ 3 




Note that (a) the Ti term is the unpolarized differential decay rate and the other three 
terms Tj [i = 2, 3, 4) are polarization-dependent; the T 2 term is due to longitudinal polar- 
ization, the 1^3 term due to transverse polarization and the r 4 due to normal polarization, 
and (b) the transverse and normal components after integrating the differential decay 
rate over the azimuthal angle G or <b v vanish, which is consistent with the so-called "null 



transverse-polarization theorem" |12| . Clearly, to utilize the polarization-dependent terms 
effectively, the initial r beam should be highly polarized, that is to say, P should be rel- 
atively large. Deferring the potential impact of the value of P to our later discussion, 
we assume P=l for the time being. In this case, each dYi/d<& 3 (i — 1 to 4) term can be 
decomposed into a CP-even part Ej and a CP-odd part Af. 

— = -(£; + A,). (21) 

The four CP-even parts Sj and four CP-odd parts Aj can be easily identified by use 
of the CP relation (|IJJ) between the r~ and r + decay helicity amplitudes and they are 
expressed in terms of the r T helicity-dependent terms dTi/d§ 3 and dTi/d§ 3 as 





d(T 1 + T 1 ) 
d$s ' 


s 2 = 


d(T 2 - r 2 ) 


s 3 = 


d(r 3 - f 3 ) 

d$s ' 


E 4 = 


d{T 4 + r 4 ) 


A 1 = 


d(Ti - f i) 


A 2 = 


d(T 2 + f 2 ) 


A 3 = 


d(r 3 + r 3 ) 


A 4 = 


d(T 4 - f 4 ) 



(22) 
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where we have used the same kinematic variables {q 2 , 0, 8} for the dTi/d&3 except for the 
replacement of by — 0. We have numerically estimated the four CP-even terms in the 
allowed phase-space points and have found that the other three Ej (i = 2, 3, 4) terms are 
negligible compared to the Si term. Therefore we neglect those terms in the following. 
The CP-even Si term and the CP-odd Aj (i — 1 to 4) can be obtained from the t t decay 
helicity amplitudes and their explicit form is listed in Appendix A. All the CP-odd terms 
are proportional to the imaginary part of the parameter £ in eq. ( |T3"D . 

An appropriate real weight function Wi(Q]q 2 ,9,4>) is usually employed to separate 
the Aj contribution and to enhance its analysis power for the CP-odd parameter Im(£) 
through the CP-odd quantity: 

(wAi) = J [wAi] d$ 3 . (23) 

of which the analysis power is determined by the parameter 

(wAi) 



(24) 



For the analysis power £j, the number Ni of the r leptons needed to observe CP violation 
at the l-o" level is 

* = shi • (25) 

where Br is the branching fraction of the relevant r decay mode. Certainly, it is desirable 
to find the optimal weight function with the largest analysis power. It is known that, 
when the CP-odd contribution to the total rate is relatively small, the optimal weight 
function is approximately given by 



<t = ^- (26) 



We adopt these optimal weight functions in the following numerical analyses with several 
concrete models beyond the SM introduced in the following. 

Although there is no CP violation in the r decays within the SM, it is possible to 
conceive several new sources of CP violation in the r decays. Among them we will consider 
models with new scalar-fermion interactions, which still preserve the SM gauge symmetries 
and have only the massless and left-handed neutrinos as in the SM. In this case, only 
four types of scalar-fermion interactions can contribute to the Cabbibo-suppressed decay 
t -> (Kir)u T [0; the MHD model |T5] and three SLQ models JTB . 



In the MHD model CP violation can arise in the charged Higgs sector with more 



than two Higgs doublets |L7j and when not all the charged scalars are degenerate. As 
in most previous phenomenological analyses, we also will assume in this MHD model 
that all but the lightest of the charged scalars effectively decouple from fermions. The 
effective Lagrangian for the decay r — > Kttu t in the assumption is then given at energies 
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considerably low compared to M H by 



C 



MHD 



2V2G F sin 9 



c 



m T m s 



X*Z(s r u l )(u tl t r ) + (—) Y*Z(s l u r )(u tl t r 



+ h.c, 



(27) 



where X, Y and Z are complex coupling constants which can be expressed in terms of 
the charged Higgs mixing matrix elements. From the effective Lagrangian, one obtain for 
the MHD CP-violation parameter Ihi(£mhd) 



Im(£ 



MHD J 



m x 



m. 



m s -m u , 



Ml 



lm(XZ*) + { — ) lm(YZ*) 
m. J 



(2f 



On the other hand, the effective Lagrangians for the three SLQ models [Tj] contribut- 
ing to the decay r — > Kirv T are written in the form after a few Fierz rearrangements: 



C 



I 


223 ^13 


SLQ — 




II 


2/232/13 r 


SLQ — 


2M| 2 [ 




, ?/23?/l*3 / 



, , 1 

[s L u R )(f c R v c TL ) + -(s L a^u R )(f c R a^ c TL ) 



sl1hU L ){t c l ^v c ) +h.c 



h.c, 



r in 

^SLQ 



2M 



(29) 
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Here the coupling constants x\j , and = 1,2,3) are in general complex so 

that CP is violated in the scalar-fermion Yukawa interaction terms. The superscript c 
in the Lagrangians £slq an d ^slq denotes charge conjugation, i.e. ip%i = il 1 2 i , Rt L m 
the chiral representation. Although the tensor parts as well as the scalar parts appear in 
Model I and Model II, we do not have the tensor contributions to r — > (Ktt)u t because we 
concentrate on the vector and scalar resonance contributions in the present work. In the 
approximation that all the CP-even contributions from new interactions are neglected, 
the size of the SLQ CP-violation effects is dictated by the CP-odd parameters 



m\ 



Im [£23^13] 



SLQ J 



(m s 



- m u )m T 4V2G F sin OqM^ 
Im [2/232/13] 



m 



im(e ; 



slq; 
slq) 



(m s - m u )m T A^/2G F sin 9 C M^ 2 

0. 



(30) 



This approximation is justified because the contributions from new physics are expected 
to be very small compared to those from the SM. 

The experimental constraints on the CP-violation parameters in (|28[) and ( pOf) depend 
on the values for u and s current quark masses, which are not well determined. Here, we 
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use for the light u and s quark masses m u = 5 MeV and m s = 320 MeV, which satisfy 
the mass relation m s — m u = 7 fx* @ • Inserting the quark mass values into (p8|) and (|30|) 
yields 



Ihi^mhd) ^ -3.2 x 10' 



2 



lm(XZ*) + — Im(YZ*) 



64 



(31) 



M&q) -39.3 Im[x 23 ^] , 



2 



Me£q) ^ " 39 - 3 (^J Im Wl, (32) 

where the W^-boson mass Mw is retained to show the Mw dependence of the parameters 
explicitly, but 80 GeV || will be used for the W^-boson mass in the actual numerical 
analysis. 

The couplings X, Y and Z in the MHD model can be constrained through the processes 
such as 5-meson semileptonic decays. Because these experimental constraints have been 
extensively reviewed in Ref . , we simply follow the analysis from which the combined 
constraint on Iul(£mhd) is obtained to be 

M£mhd)| <°- 48 > ( 33 ) 

when Mh is set to be 45 GeV. Although there are at present no direct constraints on the 
SLQ CP-odd parameters in (ISO), a rough constraint to the parameters can be provided 



on the assumption |L8| that |x' 13 | ~ \xis\ and \y' 13 \ ~ [2/13 1 , that is to say, the leptoquark 
couplings to quarks and leptons belonging to the same generation are of a similar size; 
then the experimental upper bound for the DD mixing yields 

M£slq)I < °- 15 > M£slq)I < °- 14 . (34) 

which are stronger than the constraint (^) on the MHD CP-odd parameter Iul(£mhd). 
Based on the constraints ( |33"D and ( |34] ) to the CP-odd parameters, we quantitatively 
estimate the number of the Cabbibo-suppressed r decays to detect CP violation for the 
maximally-allowed values of the CP-odd parameters: 

Im(£MHD) = 0.48, Im(£| LQ ) = 0.15, Im(£sL Q ) = 0-14. (35) 

As shown in eq. ( p5j) the branching fraction for the Cabbibo-suppressed r decay needs 
to be known beforehand, and experimentally, it is crucial to reconstruct the r momentum 
direction. So we consider the most easily-identifiable sequential decay channel r — > (K° — ► 
7t~7t + )tt~v t where an efficient identification of the three charged-pion vertex and the t + t~ 
production vertex can be used to determine the momentum of the mother r lepton. The 
branching fraction for the clean sequential decay mode is approximately 0.33% ||. Table 1 
shows the numbers iVj of r leptons required to detect CP violation through the CP-odd 
rate and polarization asymmetries A, at the l-er level for the maximally- allowed CP-odd 
MHD and SLQ parameters fl35|) and Br = 0.33%. All the CP-odd rate and polarization 
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asymmetries require a similar number of r decay events. But, we note that the results 
for the CP-odd polarization asymmetries have been obtained for completely-polarized r 
leptons. Therefore, in a realistic experiment with partially-polarized r leptons the analysis 
power of the polarization-dependent observables will be reduced. In his recent works [[|, 
Tsai has claimed that r polarization can play a crucial role in probing P, CP and T 
violation in r decays. However, we see that at least in the Cabbibo-suppressed r decay 
t — > (Ktt)u t it is crucial to highly polarize the r leptons to fully utilize the CP-odd 
polarization-dependent observables. Numerically, we find that CP violation in the MHD 
model and the two SLQ models may be detected with about 10 6 and 10 7 r leptons for 
the maximally-allowed CP-odd MHD and SLQ parameters, respectively. 



Table 1. The number of t leptons, N i; needed for detection with the at the la level, are 
determined for fx* = 45 MeV with Ihi(£mhd) = 0.48 in the MHD model, and Im(^g L q) = 
0.15 and Ihi(£|lq) — 0.14 in the two SLQ models. 



Model 


iVi 


N 2 


N 3 


iV 4 


MHD 
SLQI 
SLQII 


3.19 x 10 5 
3.26 x 10 6 
3.75 x 10 6 


3.20 x 10 5 
3.28 x 10 6 
3.76 x 10 6 


3.19 x 10 5 
3.27 x 10 6 
3.75 x 10 6 


3.21 x 10 5 
3.29 x 10 6 
3.78 x 10 6 



In summary, we have investigated CP violation from the MHD model and SLQ models 
in the Cabbibo-suppressed decay r — > Ki\v T with polarized r leptons to which both the 
J p = 1~ resonance K* and the J p = + resonance Kq contribute. In addition to the 
CP-odd rate asymmetry, r polarization enables us to construct three additional CP-odd 
polarization asymmetries whose magnitudes depend crucially on the Kq decay constant, 
fx*- Taking a QCD sum rule estimate of fx* = 45 MeV and the present experimental 
constraints on the CP-odd parameters into account, we have quantitatively estimated 
the maximally- allowed values for the CP-odd rate and polarization asymmetries in the 
MHD and SLQ models consistent with the SM gauge symmetry with massless left-handed 
neutrinos. We have found that the CP-odd rate and polarization asymmetries are of a 
similar size for highly-polarized r leptons and, for their maximally- allowed values, new 
scalar-fermion interactions may be detected with about 10 6 or 10 7 r leptons at the la 
level. Consequently, we conclude that since the r leptons of the order of 10 7 or more are 
expected to be produced yearly at the planned B factories JTT|] and the proposed r-charm 
factories ED | , it is important to look for CP violation in the Cabbibo-suppressed r decay 
t — > KttUt- even without polarized r leptons. 
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Appendix 

A. CP-even and CP-odd observables 

In the appendix the explicit form of the CP-even term Ei and the CP-odd terms Aj (i — 1 
to 4) is presented: 



E = F(q 2 ) 


2C 2 K m 2 r 











\F S \ 2 + 8P 2 K \F K \ 2 1 + 



mz 



1 cos 2 d 



l 2 f m K* \ 

-S( > K -=p f-^- + Re(0 ] Re(F K F* s ) costf 



Al = -8F{q 2 )C K m T P K 



costi1m(£)Im(F K F£) , 



(36) 
(37) 



A 2 = -8F(q 2 )C K m T P K lm(0 



cos 9 + 



m 



^ - l^J cos0costf jlm(F K i^) 



+ sin 6 sin 9 sin <pRe(F K Fg) 



(38) 



A 3 = -8F(q 2 )C K m T P K lm(0 



sin 9 cos 0+ ^"~Ff ~~ 1^ sin© cos $j lm(F K Fg) 



— cos sin 9 sin (pRe^FxFg) 

A 4 = ~8F(q 2 )C K m T P K \m(i) sin9 sm<p\m(F K F* s ) 

+ (cos sin 9 cos — sin cos 9)Re(F K Fg) 

where the overall function F(q 2 ) is given by 

G 2 F m T sin 2 9, 



(39) 



(40) 



F(q 2 



q 



2 7 7T 4 



m 



(41) 
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